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ABSTRACT 

In this work we develop a new method to turn a state-of-the-art hydrodynamical cosmolog- 
ical simulation of galaxy formation (HYD) into a simple semi-analytic model (SAM). This 
is achieved by summarizing the efficiencies of accretion, cooling, star formation, and feed- 
back given by the HYD, as functions of the halo mass and redshift. The SAM then uses these 
functions to evolve galaxies within merger- trees that are extracted from the same HYD. Sur- 
prisingly, by turning the HYD into a SAM, we conserve the mass of individual galaxies, with 
deviations at the level of 0.1 dex, on an object-by-object basis, with no significant system- 
atics. This is true for all redshifts, and for the mass of stars and gas components, although 
the agreement reaches 0.2 dex for satellite galaxies at low redshift. We show that the same 
level of accuracy is obtained even in case the SAM uses only one phase of gas within each 
galaxy. Moreover, we demonstrate that the formation history of one massive galaxy provides 
sufficient information for the SAM to reproduce the population of galaxies within the entire 
cosmological box. The reasons for the small scatter between the HYD and SAM galaxies are: 
a) The efficiencies are matched as functions of the halo mass and redshift, meaning that the 
evolution within merger- trees agrees on average, b) For a given galaxy, efficiencies fluctuate 
around the mean value on time scales of 0.2-2 Gyr. c) The various mass components of galax- 
ies are obtained by integrating the efficiencies over time, averaging out these fluctuations. We 
compare the efficiencies found here to standard SAM recipes and find that they often deviate 
significantly. For example, here the HYD shows smooth accretion that is less effective for low 
mass haloes, and is always composed of hot or dilute gas; cooling is less effective at high red- 
shift; and star formation changes only mildly with cosmic time. The method developed here 
can be applied in general to any HYD, and can thus serve as a common language for both 
HYDs and SAMs. 

Key words: galaxies: evolution - galaxies: formation - galaxies: haloes - large-scale structure 
of Universe. 



1 INTRODUCTION 

The formation and evolution of galaxies within our Universe is a 
complicated process that combines two very different mechanisms. 
On the one hand, the hierarchical growth of dark-matter structure 
drives the aggrega tion of galaxies, on time- scales that are propor - 
tional to redshift dPress & Schechtejll974l : lLacev & Cold 1 19931) . 
On the other hand, the baryonic physics determines the interplay 
between gas and stars, on time scales that are affecte d by local pro- 
cesses of cooling, star formation (SF) and feedback dWhite & Rees l 
1978l : lDekel & Silkl 19861 : [white & Frenkll99ll : lMcKee & Ostrikej 
2007h . The combination of these two disciplines shapes the com- 
plex evolution of galaxies over cosmic time. 

Models that take into account the above processes differ in 
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their level of complexity, and in the typical scales that are be- 
ing resolved or properly modeled. In general, a simple distinc- 
tion can be made between two different approaches, namely hydro- 
dynamical simulations (hereafter HYDs), and semi- analytic mod- 
els (SAMs). HYDs try to follow the evolution of a galaxy, by 
modelling in great detail the hydrodynamics and gravitation laws 
that are in play. These models often use more than 10 6 parti- 
cles to describe one galaxy, and thus allow its detailed struc- 
ture to be explored. However, HYDs are still limited to a fi- 
nite resolution, which does not allow all the processes men- 
tioned above to be followed properly. Consequently, HYDs rely 
on 'sub-grid' analytical laws, that describe SF, feedback, and 
the structure of the inter- stel lar medium (ISM ) . For a few ex- 



amnles of HYD studies, see Katz et al. ( 1996b: Governato et alj 


(19991): SDrinsel & Hernauist d2003l): 


Scannapieco et al.1 d2009l): 


ISchave et all J2010h: lAsertz et all (1201 1 


). 
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A different approach, adopted by S AMs, is to treat each galaxy 
as one unresolved object, using integrated properties to describe the 
mass of stars, cold gas, hot gas, and the black hole. Since each com- 
ponent of the galaxy is represented by one number, the dynamics 
within the galaxy is not resolved, and one needs to come up with 
laws for star formation, cooling, and feedback that are valid on av- 
erage for the entire galaxyQ Due to their simplicity, S AMs can pro- 
vide a statistical sample of galaxies, and can explore a large portion 
of their parameter space. F or more details , the reader is referred to 
some recent SAM studies: [Monaco et al. (2007); Somerville et al, 



20081 ) : iBenson & Bowed d2010h : iGuo et all (1201 ll) : IWang et al 
201 ll) : lKhochfar et all ( 1201 ll) 



In the last two decades, HYDs and SAMs have been used as 
the two major tools for studying the formation and evolution of 
galaxies. Detailed comparisons between the two approaches are 
thus important both for developing better models, and for hav- 
ing a common language to interpret different models. Following 
this reasoning, various comparisons between the two methods were 
made to date. Most of these studies have focused on the processes 
of accretion and cooling, finding some agreement at low redshift, 



and la rg er deviations a t high- z. For more d e tails, see Benson et al 



d20Qlh : lYoshida et ail d2002h; iHellv et all (120031): ICattaneo et al 
(120071) : Iviolaetall d2008l) : ISaroetal.1 b010h : LuetalJ (1201 ll) 



iHirschmann etal.1 ( 201ll) . In each of the above works, both the 
SAM and the HYD are adopting a specific model with a given 
parametrization. Thus, it is not clear if the discrepancies found 
between the HYD and SAM galaxies are due to the limitation of 
each methodology, or are just a simple outcome of the specific 
model chosen. A few other works hav e tried to quantify the physics 
of HYDs without using a SAM (e.g. iHernauist & Springel| !2003 



iRasera & Tevssie rll20o3 : iDave et al.ll201ll) . Although such studies 
are important for understanding the physics of galaxy formation, it 
is difficult to estimate the accuracy of these models for individual 
objects. 

Recently, IStringer et al.l d2010h have tried a different path to 
attack this issue, by trying to tune a SAM according to the physics 
of a HYD. These authors managed to modify a SAM based on 
iBower et"aD (120061) . so that it will roughly reproduce the history of 
one disk galaxy within a HYD. Since their work was based on only 
one galaxy, and since some deviations between their SAM and the 
HYD remained, it is still not clear how well the two methodologies 
agree. 

In this work we would like to take this approach one step fur- 
ther. We will develop a method to extract the physics of a HYD us- 
ing the simulation output, in a way that can be used within a SAM. 
We use a large cosmological hydrodynamical simulat ion, based on 
state-o f-the-art physical modelling, as developed by ISchave et al.l 
(l2010l) . Our task is to explore the level of complexity needed by a 
SAM in order to follow accurately the formation histories of galax- 
ies as modeled by the HYD within a large cosmological box. 

This paper is organized as follows. In section[2]we describe the 
HYD and the SAM used here, and the method being used to extract 
SAM ingredients out of the HYD. These ingredients are presented 
in section [3] emphasizing the differences in comparison to stan- 
dard SAMs. The galaxies produced by both models are compared 
in section HI A model with one gas phase is presented in section [5] 
showing a similar match to the HYD as in the case of the standard 



1 In more detail ed SAMs, that model e.g. the SF rate as a function of the 
disk radius I Dutton & van den Bosch 2009; Fu et al. 2010) one needs to as- 
sume an ad-hoc density profile within the disk. 



SAM. In section [6l we further discuss a few additional tests of the 
formalism, and try to pin down the reasons for its success. Lastly, 
we summarize and discuss the results in section [7] 



2 METHODS 

2.1 The hydrodynamical simulation (HYD) 

In this work we use a cosmological hydrodynamical simulation 
(HYD) based on the Overw helmingly Large Simulations (OWLS) 
project dSchave et al.ll20ld) . This project includes a large set of 
HYDs with various di fferent physical ing r edients that we r e stud - 
ied extensively by e.g. ISales et al. J2OI0I); IWiersma et al] (l201ll) : 



Ivan de Voort et all (1201 ll) : iMcCarthv et al.l ( 1201 ll) . Herewe only 
use one si mulation setup , the s ame as the 'reference model' de- 
veloped by ISchave et al. (|2010|). In brief, this model includes ra- 
diative cooling based on lwiersma et"aD (120091) . following the con- 
tributions from 11 different elements that are released by stellar 
winds from massive st ars, AGB stars and by supernovae of types la 
and II, as described in lWiersma et al] d2009h . The SF law is guided 
by the observed Kennicutt- Schmi dt law, implemented in the form 
of a pressure law as described in ISchave & Dalla Vecchial d2008l) . 
Supernova (SN) feed back is modeled by injecting SN energy in ki- 
netic form, following iDalla Vecchia & Schavd (120081) . This model 
includes neither active galactic nuclei (AGN) nor black holes. 

We ran a new simulation that is identical to the OWLS refer- 
ence model with a box size of 100 /i -1 Mpc, and 2 x 512 3 particles 
of dark-matter, gas and stars. The simulation outputs were saved 
in 68 snapshots, more than the original run, and approximately 
spaced by 200 Myr, from z — 20 to z = 0. The dark-matter particle 
mass equals 4.06 x 10 8 /i _1 M©, and baryonic particles have initial 
mass of 8.66 x 10 7 /i _1 M©. The comoving (Plummer-equivalent) 
gravitational softening is 7.8 /i -1 kpc (with a maximum value of 2 
h~ 1 kpc in physical units). The underlying cosmological parame- 
ters are: (Q m , a 8 , n s , Q b , h)= (0.238, 0.74, 0.951, 0.0418, 0.73), 
consistent with the WMAP 5-year data dKomatsu et al. 112009b . 

On each output snapsh ot we have run the FOF algorithm 
with a linking length of 0.2 dDavis etal.lll985l) to identify haloes 
with more than 20 da rk-matter particles. The SUB FIND algorithm 
dSpringel et all 1200 ll) was then used to identify subhaloes with 
more than 20 particles within haloes (i.e. the minimum subhalo 
mass ranges between lxlO 9 to 8xl0 9 /i _1 M©, depending on 
which particles are included). In our imple mentation, SUB FIND 
uses both dark matter and baryonic particles dDolag etal.ll2009l) . 
Since satellite subhaloes within a dense environment are often be- 
ing stripped of their dark-matter, subhaloes occasionally host only 
star and gas particles. In addition, fragmentation might happen 
within haloes, creating new subhaloes, with only gas and star par- 
ticles. We have con structed merger-trees of subhaloes in the same 
way as described in lSpringel et al] d200ll) . The trees include infor- 
mation on the subhaloes and their host FOF groups. 



2.2 The semi-analytic model (SAM) 

Here we describe the semi-analytic model (SAM) used in this work. 
For more details on the model, including variou s speci fic scenarios 
for galaxy evolution, see Neistein & Weinmannl d2010l) . The model 
follows galaxies inside the complex structure of merger-trees, and 
uses simple laws for cooling, SF, accretion, merging, and feedback. 
Unlike other SAMs, these laws are simplified to be functions of 
only the host subhalo mass and redshift. 
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2. 2. 1 Quiescent evolution 

Galaxies that do not experience merger events are termed to evolve 
'quiescently'. Each galaxy is modeled by three phases of baryons, 



a galaxy : (m sta r, m co id, m ho t) • 



(1) 



The definitions of m co id and mhot are motivated by the HYD: 
racoid is the mass of the cold and dense gas that is able to form 
stars (temperature smaller 10 5 K, density larger than 0.1 cm 3 ), mhot 
is all the rest of the gas within the host subhalo, including gas that 
was previously inside the subhalo but was later ejected. The exact 
definitions of the different gas phases are given in the next section. 
We note that there exist various different definitions for cold and 
hot gas in the literature. Although our definition agrees with the ap- 
proach adopted by SAMs, it is different from recent studies based 
on HYDs, as will be discussed below. In addition to our standard 
model, we will test various scenarios with differing number of gas 
phases. 

In the following we lay out the basic set of differential equa- 
tions that describe the evolution of these phases using a small set 
of a priori physical assumptions. These equations have been the 
basis for t he standard paradigm of galaxy formation for oyer 30 
years now feees & O striked! 19771 : Isilkll 19771 : 1 White & Reel 19781 : 
IWhite&Freriklll99ll) . 

A fresh supply of gas into the galaxy is provided by smooth 
accretion along with the growth of the host dark-matter subhalo. 
The efficiency of hot accreted gas is modeled by 



[m ho t J 



accretion 



U-M h ifM h >0 
otherwise 



(2) 



Here Mh is the subhalo mass, defined as the total mass of dark- 
matter particles within the subhalo. Mh is the growth rate of dark- 
matter coming from particles that are not included in other sub- 
haloes (not mergers). Square brackets are used to identify individ- 
ual processes, in this case it is the contribution to mhot due to ac- 
cretion. We allow / a to be a fun ction of the halo ma ss and redshift, 
although in standard SAMs fe.g. lCroton et al.ll 2006) it is a constant 
that equals the universal baryonic fractionQ / Ct rn = 0.1756. In 
general, a similar component of cold accretion might exist. How- 
ever, as will be discussed below, cold accretion is negligible due to 
our definition of cold gas, which includes a threshold in density. 
Hot gas may radiate and cool according to 



fc • m ho t • 



(3) 



[m co idJ cooling — L m hotJ cooling 

The cooling efficiency, f c = f c (Mh,z), is a function of the host 
halo mass Mh and the redshift only, and is written in units of 
Gyr -1 . We assume that the SF rate is proportional to the amount of 
cold gas, 

[m s tar] SF — - [m co ld] SF = fa " ^cold , (4) 

where / s = f 8 (Mh,z) has units of Gyr -1 . Gas can be heated due 
to SN explosions and move from the cold phase into the hot. In 
the HYD used here, core collapse SN events follow star formation 
after a short delay of 30 Myr. Therefore, the feedback should be 
proportional to the SF rate, 



m ho t 



: [ m coidJ feedback — 

fd [m s tar] SF = fdfs • m co ld • 



(5) 



2 For low mass haloes, reionizat ion introduces a filtering mass scale that 
gives lower baryon fractions (e.g.[spmerville 20021). 



Similar to the other ingredients, feedback is modeled by a function 
of the halo mass and redshift fd = fd(Mh,z). 

All the processes above can be united into a set of differential 
equations, 

m st ar — f s • m co id 

m co id = - (f a + fdfs) • m co id + fc • m ho t (6) 

mhot = fdfs • m co id - fc • mhot + /a ■ M h . 

Each physical process is described by one function (f x ), resulting 
in a set of linear inhomogeneous differential equations. The hot ac- 
cretion, f a • Mh, is the 'source term' that governs the total baryonic 
mass within each galaxy. The other three efficiencies (/ s , / c , fd) 
define the complex evolution of gas and stars within a galaxy. 

2.2.2 Satellite galaxies 

In this work we assume that each subhalo includes only one galaxy. 
Since subhaloes might contain only star particles, small subhaloes 
inside massive FOF groups can survive longer than in dark-matter 
only simulations. We note that although our SAM uses only the 
dark-matter mass for each subhalo, the location of the subhalo and 
its merging time are affected by the dynamical processes within 
the HYD, including contributions from the gas and star particles. 
Satellite subhaloes are defined as all subhaloes inside a FOF group 
except for the central (most massive) subhalo. Because galaxies and 
subhaloes have a one-to-one correspondence, we use the same ter- 
minology for central and satellite galaxies. 

While satellite galaxies move within their FOF group, they suf- 
fer from mass loss due to tidal stripping. This is modeled by addi- 
tional terms in the differential equations above: 

m s tar = fs ' m co ld 

m co id = - (fa + fdfs + OL c ) • m co id + fc • m ho t (7) 

m ho t = fdfs • m co id - (f c + olh) • ?™hot + /a • M h . 

The additional terms including a^, ol c are computed only for satel- 
lite galaxies, and describe the stripping of hot and cold gas respec- 
tively. In general, a similar parameter for stellar stripping can be 
added, but it is negligible in the analysis done here. Our model al- 
lows for the stripped mass to be added to the central object, or to 
be lost to the inter-galactic medium. For satellite galaxies, all the 
efficiency values f c , fd, fs are based on the subhalo mass at the 
last time the subhalo was central within its FOF group. 



2.2.3 Mergers 

In case a subhalo merges into a more massive one, we merge the 
corresponding galaxies as well, and at the same time. Mergers can 
trigger SF bursts, with an efficiency that depends on the mass ratio 
of the two galaxies: 

Am st ar = 0.56(m 2 /mi) 0-7 x m co id , (8) 

where mi , m<2 are the baryonic mass of the central and satellite 
galaxy respectively, and m co id is the sum of the cold gas masses 
of the two galaxies. T his recipe follows the resu lts of hydrody- 
namical simulations by iMihos & Hernqu ist (1994 ) and ICox et al.l 
d2008h. and was adopted by various SAMs dSomerville et alj200lF 
iCroton et al. 2006; Khochfar & Silk 2009; Neistein & Weinmannl 
l201Gh . However, as will be explained below, we do not find a strong 
evidence that these bursts are necessary to reproduce the HYD 
galaxies, and we therefore do not include bursts in our final im- 
plementation of the model. 
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2.3 How to turn a HYD into a SAM 

We would like to extract the effective laws that govern the evolu- 
tion of galaxies within the HYD. In the language of our SAM, we 
need to identify the functions f c , / s , fd , fa that summarize the pro- 
cesses of cooling, SF, feedback, and accretion. For satellite galax- 
ies, we need to determine the constants a c , ah that describe the 
stripping rates of cold and hot gas. In order to do so we follow 
each subhalo within the HYD along with its merger-tree over time, 
and keep track of all the particle information. As in the SAM, we 
assume each subhalo includes exactly one galaxy. 

At each snapshot we define three different mass components 
for each galaxy: 

• The stellar mass, ra s tar, defined as the total mass of all star 
particles within the subhalo. 

• The mass of cold gas, m co id, is the mass of all particles that are 
able to form stars within the subhalo. According to the SF law being 
used by the HYD, these are all particles with local gas densities 
larger than 0.1 cm 3 and temperatures lower than 10 5 K. 

• The component of hot gas, rahot, includes all gas particles 
that do not belong to m co id, as well as all particles that were once 
within the subhalo, but were ejected later. These ejected particles 
are assigned to the same subhalo only if they did not become part 
of other subhaloes. Note that usually in SAMs the ejected gas is 
treated as a different gas component. 

We keep track of all particles that belong to subhaloes within 
the HYD, and check which of them have changed their phase (i.e. 
^coid, rastar, rnhot) between two subsequent snapshots, or were 
accreted/stripped. For each galaxy i we define all possible transi- 
tion rates of the kind iCid^star* ^coid^hot, #hot->coid» we also 
checked that other rates, like i^hot^star are negligible. For exam- 
ple, in order to compute the SF rate we use the following sum: 



-*^cc 



h^ 1 



(9) 



Here rrij is the mass of the particle j, and the sum goes over all 
particles that started as m co id at the beginning of the time- step, 
and ended as stellar particles. At is the time in Gyr between the 
two snapshots considered. In order to compute cooling (or heating) 
rates we use a similar sum, taking into account all particles that 
started as hot (coldVparticles at the beginning of the time- step, and 
ended as cold (hot)jj 

For the accretion rate, R l _^ hot , we use the sum over all par- 
ticles that joined the hot component of the subhalo, and were not 
identified inside other subhaloes before. In addition, we take into 
account the mass of particles that were exchanged between sub- 
haloes that belong to different FOF groups. This means that parti- 
cles that are stripped into a different FOF group are subtracted from 
the accretion rate. On the other hand, particles that join the central 
subhalo coming from satellite subhaloes within the same FOF group 
are not accounted for in accretion rates. 

Whenever we have a merger event, we first sum up the com- 
ponents of the progenitor galaxies, and only then compute the rates 



3 Multiple transitions of the type hot— ^cold— ^hot might exist between two 
snapshots, but these are negligible according to the cooling time-scales that 
will be shown below. In case multiple transitions exist, not including them 
will modify the rates we measure. However, these modifications should not 
change the mass components of the SAM galaxies. This issue reflects the 
inherent degeneracy of the model equations. 



for the remnant galaxy. This means that our rates reflect the qui- 
escent evolution only, and do not include mergers explicitly. How- 
ever, mergers might still induce bursts both in the HYD and SAM, 
following, e.g., Eq.[8] Mergers can also affect other processes indi- 
rectly, like heating, cooling or SF within the HYD. This issue will 
be discussed below. 

The efficiencies for each galaxy i are defined by normalizing 
the rates: 



fa 



RU 



mi 



RL 



R 



cold— >-stars 



/a 



R. 



(10) 



(11) 



(12) 



(13) 



cold— >-stars 



In order to obtain the global efficiency law, for the full cosmological 
box, we consider only central subhaloes within their FOF groups. 
We then split the sample of subhaloes into bins of different mass 
and redshift. For each bin the average efficiency is defined by aver- 
aging the nominator and denominator separately. For example, 



fs(M h ,z) 



r) 



(™cold) 



(14) 



Here averaging is done over all galaxies within the same Mh and z 
bin. Quite arbitrarily, we choose bins of 0.2 dex in Mh, and 7 bins 
in cosmic time, spaced by ~ 2 Gyr. We have checked that finer 
bins do not modify the results of this work0 The bins in cosmic 
time are much wider than the time between two subsequent snap- 
shots. Consequently, the time average typically includes 10 differ- 
ent snapshots. 

When computing the accretion efficiency, we use the fact that 
within the SAM, negative dark-matter accretion events are treated 
as zero, and are not inducing negative gas accretion. To make this 
approach consistent with the average value of f a measured from 
the HYD, we set all negative values of Mh to zero first, only then 
do we average Mh and compute f a : 



f,{M hl z) 



(-f^-s-hot) 



(max (Mh,o)) 



(15) 



This way of averaging guarantees that the total baryonic mass 
within our SAM galaxies will agree with the HYD. 

We have saved stripping and accretion rates for satellite galax- 
ies, and recorded the amount of mass flowing into the central sub- 
halo, in comparison to the total mass being stripped. In general, the 
stripped mass is best described by a normalized efficiency, i.e. the 
ratio ah = Arahot/^hot, as was defined in Eq. [7] However, we 
found that both ah and a c are changing as a function of the host 
subhalo mass and redshift. In addition, the efficiencies of feedback 
and cooling for satellite galaxies are somewhat different than those 
for central galaxies. This partially depends on the definitions of the 
various gas phases, and how each is being stripped. In this work we 
have chosen the simplest model possible, taking into account only 
constants a c and ah- Furthermore, we use the same cooling and 
feedback efficiencies as for central galaxies. We will show below 



4 The SAM used here (Neistein & Weinmann 2010) automatically inter- 
polates the input values of / s , fa, fc, fd into a fine grid in halo mass and 
time. 
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Lookback time [Gyr] 



Figure 1. The efficiencies of cooling, accretion, SF, and feedback following 
the main-progenitor history of a specific galaxy within the HYD, hosted by 
a subhalo of mass ~ 10 14 Mq at z = 0. Solid lines show the different 
efficiencies for this galaxy as measured from the HYD at all snapshots. 
The thick dashed lines correspond to the average efficiency for all galaxies 
within the HYD, including only central galaxies with the same host subhalo 
mass and at the same redshift as the galaxy plotted in solid lines. The thin 
dashed lines show the standard deviation for the same sample of galaxies 
above. The mass components of this galaxy are shown in Fig. [2 



that this solution is reasonably accurate for satellite galaxies. We 
plan to investigate this issue more closely in a future work. 

In the hydrodynamical simulation, star particles can lose some 
of their mass due to stellar winds and SN. This mass loss is com- 
puted using a stellar populati on synthesis model and is added to 
the surrounding particles fsee lWiersma et al.l Eo09). Although this 
process can be easily modeled within the SAM, it complicates the 
interpretation of the results. This is mainly because the stellar mass 
loss at a given epoch is the outcome of the SF history over a few 
Gyr. Therefore, the rates measured from the HYD would not be 
instantaneous, and might include less scatter with respect to the 
SAM. We therefore assume that all particles have a fixed mass, and 
compute all rates and efficiencies using this assumption. This as- 
sumption is also being used when comparing the results of the HYD 
against the SAM. Consequently, the total baryonic mass within 
galaxies is sometimes higher than the universal baryonic fraction. 
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Figure 2. The baryonic components for the main-progenitor history of a 
specific galaxy. Solid lines show the mass of stars, cold gas, and hot gas 
as measured from the HYD. Dashed lines represent the same components 
within our SAM, using the same dark-matter subhaloes (see section [4} . The 
efficiencies of cooling, accretion, SF, and feedback for the same galaxy are 
plotted in Fig.[T] 



For completeness, we show in the Appendix all the efficiencies 
from the HYD when using the proper mass for each particle, as 
was used in the simulation. 

In Fig. Q] we show the different efficiencies for the main- 
progenitor histor)|3of one massive galaxy within the HYD, in com- 
parison to the global averages using all central subhaloes of the 
same mass and time within the HYD. It seems that the random- 
ness in the efficiencies of one galaxy is not too big, and is averaged 
out over time (except for a few narrow peaks that should not affect 
the masses of stars and gas significantly). For example, f c and f a 
show deviations on time- scales of one snapshot (200 Myr), with 
no significant trends over larger time- scales. On the other hand, fa 
and f s show deviations from the average efficiencies that are last- 
ing for ro 2 Gyr. Overall, the behaviour of one galaxy seems to be 
very regular, and does not show significant deviations larger than 
the standard deviation (STD) computed using all the galaxies in the 
HYD. The total mass in stars, cold gas, and hot gas for the same 
galaxy are plotted in Fig. |2] We will show below that once we use 
the SAM over the same merger-trees, the agreement between HYD 
and SAM is very good, also when comparing individual objects. 



3 THE PHYSICS OF THE HYD 

The different efficiencies extracted from the HYD should describe 
the various physical processes involved in forming galaxies. As we 
will see below, they allow the SAM to accurately reproduce the 
population of galaxies in the HYD. This means that we have a reli- 
able estimate of the net effect that heating, cooling, accretion, and 
SF have on galaxies within the HYD. 



5 The main-progenitor history is defined by following back in time the 
most massive progenitor in each merger event. Note that at high redshift, 
the subhalo that belongs to the main-progenitor branch might not be the 
most massive within its merger tree. 
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Figure 3. The smooth accretion rate of baryons as derived from the hydro- 
dynamical simulation. / a is defined as the ratio of the smoothed hot accre- 
tion, over the dark-matter smooth accretion, averaged in bins of halo mass 
and time. Each solid line represents a different redshift bin, all the other 
lines are shown only for z = 1: The dashed line shows the average plus one 
standard deviation in / a ; The dotted-dashed line shows the average minus 
one standard deviation, after averaging out / a for all the progenitors within 
each tree, and at ~ 10 different snapshots (all the merger- trees are rooted 
at z = 0); The dots represent the average plus one standard deviation of 
f a after averaging over different progenitors within a tree, but not within 
different snapshots. The thick dashed line is the universal baryonic fraction, 
VL b /(VLm -Q b ) = 0.213. 



3.1 Smooth accretion 

The values of the accretion rate, / a , that were extracted from the 
HYD are shown in Fig. [3] We plot only the 'hot accretion' com- 
ponent as we do not detect an accretion of cold gas into galaxies. 
Alth ough this seems to be in conflict with various recent studies 
(e.g.lKeres etalJl200l 120091 : iDekel et alJEooi Ivan de Voort et all 
2011) it is a result of the different definitions of 'cold gas' that are 
being used in the literature. Here we define cold gas as the gas that 
is able to form stars, requiring high densities (larger than 0.1 cm 3 ), 
and not only low temperatures. This is a different definition from 
most other studies based on HYDs, that often define a gas particle 
to be cold if it was not previously heated to the virial temperature of 
its halo. Here we adopt a more straight-forward definition of cold 
gas, based on the SF law. Using our definition, there is no evidence 
for 'cold accretion' at all redshifts and for all subhalo masses. This 
fact is reasonable, because star-forming gas might form stars before 
it joins the subhalo, and will therefore be identified as a separate 
galaxy. 

Unlike in standard SAMs, where / a is assumed to be a con- 
stant, here f a shows a significant dependence on the subhalo mass, 
decreasing by a factor of ~ 10 from subhalo masses of 10 13 to 
1O 11 M at z — 0. This is surprising, considering the fact that all 
galaxies used here are the central objects inside their FOF groups. 
We have checked this effect further, and tested a model in which 
negative gas accretion (stripping) is allowed whenever the dark- 
matter mass decreases. Using this new assumption, the accretion 
efficiencies become much clo ser to a constant, with dev iations of a 
factor of ~ 2, consistent with lvan de Voort et al.l J2Q1 lh (the equiv- 
alent plot of f a in this case is shown in the Appendix, Fig. lAlb . Our 
conclusion is that the low accretion rates shown here for low-mass 



Figure 4. The cooling efficiencies extracted from the HYD, and averaged 
as a function of halo mass and time. f c is defined as the ratio of the cooling 
rate over the mass of hot gas. Plotted lines are using the same definitions as 
in Fig [3] Dashed lines on the bottom right are proportional to one over the 
cosmic time. 

subhaloes are compensating for the negative accretion events. Mod- 
els based on average dark-matter accretion rates should therefore 
use the efficiencies quoted in the Appendix. We have also checked 
that the total baryonic fraction within galaxies agrees with the ac- 
cretion rates given in the Appendix. In terms of the comparison 
made here, once we include a mechanism for gas stripping within 
central galaxies, following the dark-matter evolution, we do not get 
a better agreement between the HYD and SAM. We therefore adopt 
the solution of positive accretion only, as it is more simple to im- 
plement. 

As was discussed in section |231 the accretion rates shown here 
are based on a fixed particle mass, without taking into account the 
mass loss due to stellar winds and SN within the HYD. Conse- 
quently, accretion rates can have values larger than the equivalent 
cosmic value, / (Q m — Qb). In the Appendix we plot all the effi- 
ciencies using the proper mass for each particle within the HYD. 
This effect changes the overall normalization of each efficiency 
slightly, but it does not change the trends with subhalo mass and 
time. 



3.2 Cooling 

The next process for which the efficiencies are required is cooling. 
In Fig. [4] we show average cooling efficiencies for all the galaxies 
within the HYD. Here we should keep in mind that the component 
of 'hot gas' includes gas particles that were ejected out of the sub- 
halo. Therefore, the cooling efficiencies are normalized by the sum 
of both ejected and hot gas, according to Eq. QT] For a different 
definition of cooling efficiencies, where hot and ejected phases are 
treated separately, we refer the reader to the Appendix. We note 
that the dip in the cooling efficiencies seen at a subhalo mass of 
~ 7 x 10 11 M is due to the combination of both hot and ejected 
phases. 

In general, cooling efficiencies are showing a roughly con- 
stant behaviour as a function of subhalo masses for subhaloes lower 
than ro 10 12 M©, and go down for more massive subhaloes. This is 
qualitatively in agreement with semi-analytic models. However, the 
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Figure 5. Star-formation efficiencies within the HYD, defined as the ratio 
of the SF rate over the mass of cold gas. For the line definitions see Fig. [3] 



Figure 6. Feedback efficiencies measured from the HYD, and averaged as 
a function of halo mass and time, fa is the ratio between the heating rate 
and the SF rate. The thick gray line is proportional to the virial velocity of 
subhaloes to the power of -3/2, other lines are defined as in Fig. [3] 



dependence of the cooling efficiencies on cosmic time is stronger 
than a simple linear dependence. Since the dynamical time within 
subhaloes is proportional to the cosmic time, cooling cannot be 
modeled only by the infall time of gas into the centre of haloes. This 
might be a result of the cooli ng process itself, and its dependence 
on the hot gas properties (e.g. lMcCarthv et al.l2008l ; IWiersma etaD 
l2009h . 

In the terminology of 'cold accretion' mode, where all the 
accreted gas is a ssumed to be falling in narrow streams (e.g. 
iDekel etalJl2009h . the process of 'cooling' describes the time it 
takes the stream to reach the central disk, and become dense enough 
to be a part of our definition of m co id. A stronger dependence of 
cooling on the cosmic time might mean t hat accretion through fila - 
ments is more relevant at high redshift (Ivan de Voort et al.ll201ll) . 
A different option is that trajectories of streams are more ra- 
dial at high-red s hift. T his last fact was already pointed out by 
IWeinmaiinetallbOllh . and is in agreement with the orbits of sub- 
haloes within cosmo logical dark-matter simulations dWetzelll201ll : 
iHopkins et aIll201Qh . 

The cooling efficiencies shown here can be compared to stan- 



White & Frenk 


( 1991). This issue was investigated by various stud- 


ies in the past ( 


Benson et al.ll200ll:lYoshida et al.l 


20021: Hellv et all 


20031: Cattaneo et al. 2007:be Lucia et all2010t 


Crain et all2010t 


Lu et al. 2011), mostly claiming some agreement between differ- 



ent SAMs and HYDs, and so me noticeable devia tions (especially 
at high redsh ift). For exampl e JCrain et al.l d2010h showed that the 
algorithm of IWhite & Frenkl dl99lh strongly overpredict cooling 
rates due to the specific entropy profiles o f gas within haloes. In 
addition, various SAMs that are based on Whit e & Frenkl dl99lh 
can show significant deviations i n cooling rates due to the detailed 
implementation of the algorithm (iDe Lucia et alfcoiOh . 

We find significant differences w hen comparing the cooling 
efficiencies here ag ainst the one used by De Lucia & Blaizotl (12007) 



and summarized in Neistem & Weinma nn (2010). For example, at 
z — 1 JPe Lucia & Blaizot d2007h predict cooling efficiency of — 1 



3.3 Star formation 

Fig. [5] shows the SF efficiencies found within the HYD. The time- 
scales for converting the cold gas into stars range from ~ 3 Gyr for 
low mass subhaloes at z = 0, to ~ 1 Gyr for massive subhaloes at 
z = 3. The low-redshift values are roughly consis t ent with the ob- 
serva tional constraints JSchiminovich et al1l2010l : ISaintonge et al.l 
201 H). However, the dependence on redshift found here is much 
smaller than what is usually assumed in SAMs, where the conver- 
sion efficiency is prop ortional to the cosmic time (see, however, 
For example, IWang et all (EoT H) showed 
the SF efficiencies as a function of subhalo mass for various mod- 
els, where in standard SAMs the difference between z — 3 and 
z = reaches an order of magnitude. 

Interestingly, the SF efficiencies show a double power-law be- 
haviour as a function of subhalo mass, where the peak efficiency 
is located at ~ 10 12 — 1O 13 M , depending on the specific red- 
shift. For high-mass subhaloes the SF is not significantly reduced. 
Consequently, the high fraction of passive galaxies within massive 
subhaloes is not related to a reduced S F efficiency, but rather to ga s 
consumption, environ mental effects dKhochfar & Ostrik er 2008), 
or AGN feedback fe.g. lCroton e7alll2006l) . 



3.4 Feedback 

The feedback efficiencies extracted from the HYD are plotted in 



-3/2 



below 



Gyr 1 at subhalo mass of 10 11 M©, roughly a factor of 10 higher 
than what is found here. 



Fig. [6] These seem to follow a power-law of the type 
subhalo mass of - 10 13 M Q , where v is the virial velocity of the 
subhalo. Above this mass, the feedback efficiency shows a modest 
upturn. The feedback efficiency represents gas that is heated from 
the cold phase into the hot component, and possibly ejected out of 
the subhalo. Within the OWLS reference model, SN feedback is 
implemented in kinetic form using a constant wind velocity of 600 
km s _1 . This causes the feedback to be come inefficient for halo 
masses greater than a few times 10 11 M Q (Dalla Vecchia "& Schavd 
120081 : ICrain et al.ll2009l : ISchave et allEoiol : lHaasl2010h . 

SAMs usually assume a power-law efficiency wi th very dif- 
ferent in dexes. For exampl e, IDe Lucia & Blaizotl J2007b assumed a 
constant: ICole etalJbOOOh have used a power of -2, following the 
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Figure 7. Comparing the masses of individual galaxies, SAM against the 
HYD. Each panel represents a different mass component as labeled, using 
only central subhaloes at z = 0. The panels show the two dimensional 
histogram of the pairs (rasAM, ^hyd)) describing the mass of the same 
objects in both models. The pixels are color-coded according to the log 
of the number of objects. The mean difference between HYD and SAM 
is lower than 0.08 dex for all mass components. STDs are 0.08 for ra s tar 
(except for ~ 10 11 M©, for which the STD goes to 0.13 dex). For rahot 
and ratotai the STD starts at ~ 0.08 dex for low mass galaxies and reaches 
~ 0.04 for massive galaxies. The STD for m co \& is around 0.2 dex. 



potential of the host halo; and lGuo et al] (1201 lh assumed a power 
of -3.5. All these are very different from what is found here. The 
feedback efficiency within a model that includes three phases of 
gas (ra co id, wihot, and ra e j c ) is plotted in the Appendix. 



4 COMPARING MODEL GALAXIES 

In this section we compare the results of the SAM with the HYD. 
The SAM uses only the physical ingredients that were described 
above, i.e., / a , / c , / s and /a, and was run using merger-trees ex- 
tracted from the same HYD. We specifically use the same values 
plotted in Figs.[3H6j using three more bins in cosmic time. As was 
explained in section 12.31 in order to keep the model simple we do 
not attempt to model satellite galaxies with full accuracy. We set the 
values of a c and ah to and 0.3 respectively, because they provide 
an effective behaviour which is similar to that of the HYD. The 
stripped gas is added to the central galaxy within the FOF group. 
Other than that, the SAM has no free parameters, and no tuning 
was done. 

The model galaxies from the HYD and our SAM are com- 
pared in Fig. [7] By matching the same subhaloes from the HYD & 
SAM, we are able to show the agreement between the models on an 
object-by-object basis. Unlike in previous studie s that showed large 
devia tions between SAM and HYD galaxies (e.g. lHirschmann etaD 
201 H), here the two models agree quite well. For central galaxies 
the STD of differences is less than 0.1 dex, for all redshifts and 
for the various galaxy components (except for ra co id, which usu- 
ally includes only a few tens of particles within a galaxy). The total 
mass in baryons, ra to tai = ra s tar + ra co id + rahot, is shown as a 
probe of the accuracy of the accretion rates, / a . Note that the mass 
of gas particles within the simulation is 8.64 x 10 7 /i _1 M , so 
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Figure 8. Comparing the mass of individual galaxies, SAM against the 
HYD. Data was derived in the same way as in Fig. [7] but for the popu- 
lation of satellite galaxies at z = 0. Mean (STD) differences between SAM 
and HYD reach 0.1 (0.2) dex for both ra s tar and ratotai* an d 0.3 (0.5) dex 
for ra hot . 



masses below 10 10 M© include less than 100 particles and suffer 
from various numerical artifacts, also within the HYD. 

We found that deviations between ratotai in the SAM and the 
HYD correlate strongly with deviations in rahot, and are the reason 
for most of the scatter found in rahot • This is a consequence of the 
fact that most of the baryonic mass is located in rahot- A similar 
(but weaker) correlation exists between deviations in ratotai and 
ra s tar- Even though the stellar mass is affected from various addi- 
tional processes that seem to be more complicated than accretion, 
the deviations in f a between the HYD and the SAM still affect 
ra s tar- This fact can also be seen in Figs. [3]|6l In these plots we 
show the STD of each efficiency at z — 1, after averaging out all 
progenitors of the same galaxy at z = (dotted-dashed lines). It 
is evident that the scatter in f a between different galaxies is signif- 
icantly higher than the scatter in other efficiencies. This might be 
a result of different merger-histories for subhaloes that live inside 
differ ent environment densities (the assembly bias effect iGao et al.l 
120051) . 

The small scatter in ra s tar between the HYD and the SAM 
(0.08 dex) is interesting in view of the larger scatter in ra co id (0.2 
dex). Although the masses of ra co id are usually below 10 10 M©, 
and are therefore not numerically reliable, these masses are re- 
sponsible for making stars, and somehow produce a small scatter 
in ra s tar. This issue will be discussed in section [63] A different 
contribution to the scatter in ra co id comes from mergers, and will 
be discussed in section |4~T1 

We have explored the larger scatter in ra s tar at values of 10 11 
by running a different HYD, with a different feedback model (Dalla 
Vecchia & Schaye, in preparation), and by increasing the resolution 
of the efficiency bins (both in mass and time). The high-resolution 
efficiencies were not making any noticeable change, but the HYD 
with a different feedback model results in a significant smaller scat- 
ter at ra s tar = 10 11 M©. It might be that the kinetic feedback pre- 
scription used by the HYD affects the hydrodynamical state of the 
gas in a way that is different than other cooling and heating chan- 
nels. These changes might complicate the simple distinction done 
here between cold and hot gas. In addition, it might be that the 
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Figure 9. Comparing the mass of individual galaxies, SAM against the 
HYD. Histograms are the same as in Fig. [7] but for the population of cen- 
tral and satellite galaxies at z = 2. For central galaxies, the mean and STD 
of differences between HYD & SAM are similar to the numbers quoted in 
Fig. [7] For satellite galaxies the agreement is roughly two times better than 
in Fig. [8] (in dex units). 



transition between effective and ineffective feedback is sensitive to 
other properties of the subhaloes other than the subhalo mass. 

A comparison for satellite galaxies at z = is shown in Fig. [8] 
Here the deviations between the SAM and the HYD are larger than 
for central galaxies, reaching 0.2 dex for ra s tar. We have tried a 
model in which stripping of satellite galaxies follows the stripping 
of dark-matter, according to IWeinmann et al.l d2010h . This model 
did not improve the match between the HYD and the SAM, prob- 
ably because satellite galaxies experience on average different effi- 
ciencies of cooling and feedback, as was discussed in section |231 
This can be seen in Fig. [8] where ra to tai for satellite galaxies be- 
haves better than ra s tar and rahot, hinting that the total amount of 
stripping is modeled properly. The physics of satellite galaxies is 
complicated, and deserves more attention than we give it in this 
work. 

In Fig. [9] we compare galaxies at z = 2, finding similar trends 
to z = as discussed above. The mass of cold gas is much higher at 
this redshift, so the agreement and deviations are clearer. For satel- 
lite galaxies the agreement is much better than at z = 0, probably 
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Figure 10. The mass functions of ra s tar, ra C oid> and rahot in the HYD 
(solid lines) and in the SAM (dashed lines). All galaxies at z = are se- 
lected. 



because these galaxies had much less time to evolve within their 
group, and deviations have not accumulated yet. 

The mass functions for ra s tar, ra co id, and rahot using both the 
HYD and the SAM are shown in Fig. [TO] Overall, the agreement 
between the two models is very good for all mass components. 
This is expected, as the agreement for individual objects is good. 
For rahot < 10 10 M the sensitivity of the HYD to resolution ef- 
fects seems to be high. This probably hints to the dependence of the 
cooling mechanisms on resolution. Fig.fTOlincludes all the galaxies 
within each model. We note that the other figures in this section 
only show galaxies that exist both in the SAM and in the HYD. 
Due to our definitions, galaxies that just emerged within the HYD 
(and do not have any progenitors) are not included in the SAM. On 
the other hand, the SAM keeps a small population of galaxies that 
do not have descendant subhaloes. These two populations are quite 
small and do not affect the mass functions. 

The SF rates for individual objects are compared in Fig. Qj] 
Unlike the integrated properties shown above, the SF rates show 
stronger deviations between the two models, with a STD of ~ 0.5 
dex. The scatter in ra co id is about 0.2 dex, meaning that the de- 
viations in the SF rate are dominated by variations in the SF effi- 
ciency. In addition, most of the population of galaxies at low red- 
shift form stars at a low rate, ~ 1 M yr _1 . This rate corresponds 
to just a few gas particles within a snapshot, increasing the scatter 
between HYDs and SAMs. As can also be seen in Fig. [5] the scat- 
ter in f s between all galaxies is large, reaching a factor of 3 for low 
mass subhaloes. However, when / s is averaged over all the progen- 
itors within a tree, and over a few snapshots, the scatter goes down 
dramatically. Apparently, the deviations in SF rates result in much 
smaller deviations for stellar masses. We will examine this issue in 
section [631 



4.1 How important are merger-induced bursts? 

Our SAM does not include any merger-induced processes, like SF 
bursts, heating or cooling. We have tested the contribution of merg- 
ers to the models in various ways. First, we have computed the 
average baryonic efficiencies after excluding galaxies that had a 
major merger event in the last 0.5 Gyr, or that have a major satellite 
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Figure 11. Comparison of SF rates for individual galaxies within the mod- 
els, SAM against the HYD. Central and satellite subhaloes are plotted in 
dots and circles respectively. The solid line shows the values where both 
SF rates agree. Top and bottom panels show results for z = and z — % 
respectively. 



galaxy at a distance smaller than 0.5 Mpc ('major mergers', and 
'major satellite' are defined to have a mass ratio larger than 0.2). 
Applying this selection criterion does not change the results of the 
efficiencies in a noticeable way. The agreement between the SAM 
and HYD galaxies does not change either. There is only a minor 
change in the average agreement in ra to tai. However, it might be 
that the number of merger events is small, and cannot affect the full 
population of galaxies within a cosmological box. 

A second test we carried out is to run our SAM with an addi- 
tional recipe for merger-induced SF bursts. We have used the stan- 
dard recipe given in Eq.[8]above. However, in terms of the compar- 
ison made here, this recipe does not change the agreement between 
the HYD and the SAM. 

Since the effect of mergers might accumulate with time, we 
would like to define a quantity that is related to the number of merg- 
ers a galaxy had in its past. Consider a galaxy at z = 0, and define 
the mass in stars, rrn (i > 1), that was accreted from each satel- 
lite galaxy i within the merger tree. All the stars that were formed 
within the main-progenitor branch are termed mo. Using ra s tar, the 
stellar mass of the galaxy today, we define: 



y- ( mi V 

frt WstaJ 



(16) 



If the stellar mass of a galaxy is built from N s equal values 
of mi, then <j m would equal 1/y/Ns. Low values of a m indicate 
on many merger events, while high values (close to unity) indicate 
that the galaxy is built from one branch only. The meaning of using 
a m can be related to the comparison we make between the HYD 
and the SAM. Assume that each progenitor galaxy within the SAM 
includes some random, normally distributed error in stellar mass 
(with respect to the HYD), that is proportional to its mass, rm. In 
this case the relative error in the sum of all masses (i.e. the error in 
the SAM prediction for ra s tar) will equal cr m . 

We have computed cr m for each galaxy within our SAM, and 
measured the correlation between cr m and the deviations between 
the SAM and HYD galaxies. We found that the mass within the 
SAM galaxies is higher than within the HYD for galaxies with 
more mergers (galaxies with lower cr m ). This effect is true for all 
mass components (ra s tar, m co \d, and rahot), but it is strongest for 
ra co id. This means that galaxies within the HYD lose some of their 
mass in each merger event, or that satellite galaxies lose some of 
the stripped mass to the inter-galactic medium. This effect is not 
modeled by our current SAM, but it should be straight forward to 
add it, once the treatment of satellite galaxies is more accurate. 

To conclude, we do not find any significant evidence for 
merger-induced SF bursts within the HYD used here. This seems 
to b e in conflict with previous simulations o f galaxy mergers 
(e.g. lMihos & Hernauisdll994l:ICox et al.ll2008h . However, as was 
suggested by Moster et al.1 ( 201lL it might be that the presence 
of hot gas in subhaloes regulates the efficiency of bursts within 
our simulation. It might also be that the mass gained in bursts 
does not contribute mu ch to the total stellar mass with in galaxies 
(Khochfar & Silk 2006), and that mergers are too rare jLotz et all 
120081 : iHopkins etalll2010h . On the other hand, we do find sig- 
nificant evidence for ma ss loss within mergers, as was pointed 
out by previous studies (Monaco et al. 2006 ; IPurcell etaD l2007l ; 
IConrov et al1l2007l : lYang et al.ll2009h . 



5 THE ONE-PHASE MODEL 

The analysis above was based on a SAM with two different gas 
phases within a galaxy, ra co id and rahot. As an alternative, this sec- 
tion describes a model with only one phase of gas, allowing us to 
study the uniqueness of the SAM equations. The one-phase model 
has been exp lored by othe r studies in the past. In one of the earliest 
SAM works, ICold dl99lh used modelling of cold gas in haloes to 
predict the galaxy luminosity function. Although SAMs are usu- 
ally based on three phase s of gas in galaxies, t he one-phase model 



was re c ently explored bvlBouche et aL d201Qh ilKrumholz & Dekel 
(1201 ll) : iKhochfar & Silkl ( 1201 ll) : iDave etaTT(l201 lh . In what fol- 
lows, we will try to emphasize the points of similarity and differ- 
ence with respect to these previous works. 



5.1 A SAM with one phase of gas 

For each galaxy we define 

ra gas = ra ho t + ra co id ■ (17) 

Thus, ra gas includes all gas particles within the subhalo, as well 
as gas particles that were ejected from the halo. In this case, the 
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equations that govern galaxy evolution are: 

777* star — fs ' ^T'gas 
rhgas — - fs • m g as + /a ' Mh • 



(18) 



Within the SAM, these equations assume that / s and / a are inde- 
pendent of the galaxy components ra sta r and eras • 

However, if 

we compare this model to the standard model with the two phases 
above (Eq.[6]), we get that / s = f s m co id / m gas . Consequently, if the 
standard model is accurate, the one-phase model cannot be treated 
as a set of ordinary differential equations as we assume below. 

The one phase model is useful because it is simple, and in 
case f s does not depend on m gas , it can be integrated to follow 
the evolution of one object, with no mergers. First, we define the 
integral, 



P(£o,£i 



exp 



7" 



fs&t 



(19) 



where the values of f s within the integrand are computed for the 
specific mass history of the subhalo. Once P is computed, the so- 
lution of the set of equations can then be written as 



{t) = m (t )P(t ,t) + 

/a(tl)Mh(tl) 



(20) 
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It should be emphasized that this solution is valid for one branch 
only, and cannot be expanded easily to include all the progenitors 
within a merger-tree. This is because both f a and / s depend on the 
mass of each progenitor. 

A model with one gas phase can be written in various ways. 
For example, we could have a model that is based only on th e 
cold gas within a gala xy, a s was done bv | Bouche et al.l (l2010h : 
Khochfar & Silk ( 2011) and iDave et all J2Q1 lh . It is evident from 
Eqs.[6]that a straightforward way to do it is to assume that the hot 
gas component equals the total baryonic mass, neglecting the 3rd 
equation, and keeping the two first equations almost unchanged. It 
should be noted that we then get a set of equations that is differ- 
ent from the above studies. Here the source term is not f a Mh, but 
fcf&Mh, where f a represent the fraction of hot gas within the sub- 
halo. Since this approach is different from our standard model, we 
do not explore it further here. 

We have also tested a model that includes three phases of gas, 
using m e j c in addition to the two phases in the standard model. The 
mass within m e j c takes into account all the gas particles that were 
ejected from the subhalo. These are part of rahot in the standard 
two-phase model. For most of the results quoted below, this model 
shows a similar agreement to the HYD, and is thus not discussed in 
detail. We present the feedback and cooling efficiencies computed 
for this model in the Appendix. 



5.2 Star formation 

In Fig. [12] we plot SF efficiencies for the one-phase model, describ- 
ing the efficiency of transforming the total amount of gas within a 
subhalo into stars. The values of / s combine the behaviour of / s , 
/ c , and fd from the standard model into a one, compact form (see 
Eq.[[8]>. 

The values of / s shown in Fig. [12] are d ifferent from 
what is usually assumed within one-phase models |Bouche et"al 



2O10l : iKrumholz & DekefeOl lI : iKhochfar & Silkll201ll : lDave et al 



201 lh . For example, the SF efficiency derived from our HYD 
changes rapidly as a function of subhalo mass, while other works 
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Figure 12. Star-formation efficiencies within the one-phase model, / s , as 
defined in Eq.[l8] The line-types are the same as in Fig. [3] Plus symbols are 
showing the average values of / s using a sample of subhaloes with mass 
lower than ~lO 12 M0atz = O, along with all their progenitors. 



assume a constant dependence on subhalo mass. As was mentioned 
above, there are other important differences between the model 
used here and previous studies. Here the gas mass corresponds to 
all the gas inside the subhalo, including both cold and hot compo- 
nents, while previous studies have used the cold gas only. In ad- 
dition, previous models were applied to average main-progenitor 
histories, while here the model is evolved through the complicated 
structure of merger- trees. 



5.3 Galaxies in the one-phase model 

In Fig. \\3\ (upper left panel) we compare the results of the one- 
phase SAM discussed above, to the original HYD galaxies. The 
SAM is based on the efficiency shown in Fig. [12] Although the 
one-phase model is significantly simpler than our standard model, 
the results of this model agree well with the HYD, at the same 
level as in our standard model (there is, however, a small systematic 
deviation of ~0. 1 dex that is seen only in low mass galaxies within 
the one-phase model). We have checked that a similar agreement 
is obtained at higher redshifts. To summarize, we find the same 
accuracy in matching the SAM against the HYD when using one, 
two or three gas phases for each galaxy (the same is true also for 
the gas components). 

The agreement found for the one-phase model proves that the 
SAM equations are not unique, and that different models can ac- 
curately reproduce the same HYD. We will discuss the reasons for 
this behaviour, and the related implications in section [Ol 



6 MORE TESTS 

6.1 Efficiencies based on one galaxy 



The efficiencies plotted in Figs. [31161 have usually a large scatter, 
when computing the variance over the full population of galaxies 
within the box. However, we also show that once efficiencies are 
first averaged over all the progenitors of each z = galaxy, the 
variance decreases significantly. This is shown as the dotted-dashed 
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Figure 13. Comparison of the SAM against the HYD, using deviations 
for individual galaxies. The histogram were derived in the same way as 
in Fig. [7] Here we plot only stellar masses at z = 0. The upper left panel 
shows the results of the one-phase model (section |5H . Other panels show 
a comparison using all galaxies within the box, using SAM efficiencies that 
are extracted from the history of only one massive galaxy at z = 0. Each 
panel is based on efficiencies extracted from a different galaxy, randomly 
selected. 



lines in the efficiency plots. As a result, it should be possible to ex- 
tract the SAM efficiencies using the history of one massive galaxy 
within the HYD. 

Fig. [13] shows results of various different SAMs. In each 
model we extracted the SAM efficiencies using only one random 
massive galaxy within the HYD, along with all its progenitors. The 
results of these models nicely reproduce the population of the HYD 
galaxies within the full box. This test proves that one can use a 
single high-resolution zoom simulation in order to extract the net 
result of the physics of a HYD. Once the efficiencies are known, 
the SAM can use a large statistics of merger-trees, based on dark- 
matter only simulations, to explore the same physics as was used in 
the high-resolution HYD. 

At stellar masses of ~ 10 11 M© the agreement between the 
SAM and HYD is less good. As was discussed in section [4J this 
might indicate on a larger variation in feedback efficiencies be- 
tween subhaloes of the same mass. 



6.2 The dependence of efficiencies on time-steps and on 
subhalo mass 

The accuracy of the efficiencies discussed here might depend on 
the number of output snapshots used to extract information from 
the HYD. For example, since the feedback delay time is 30 Myr 
after any SF episode, and wind particles are not allowed to form 
stars for an additional 15 Myr, the snapshot spacing should not go 
below ~ 50 Myr. To explore the sensitivity of our method to the 
time between snapshots, we have tested the same method using half 
of the snapshots in our main run. The full analysis described above 
was repeated using half of the snapshots, including merger-trees, 
computation of efficiencies, and running the SAM. We did not find 
any significant changes in the results, and we conclude that our 
method is not sensitive to the specific choice of snapshot spacing. 
It might be that HYDs which include different physical ingredients 
will show more sensitivity in this respect. 



In Fig. Q2] we plot in symbols the efficiency, / s , using a range 
of low mass subhaloes (< 10 12 Mq) within the HYD. We select 
subhaloes at z — and use all their progenitors for computing the 
average / s . It is clear that the average / s for these subhaloes agrees 
well with the total average within the box. This point allows us to 
use our method in order to increase the dynamical range of HYDs. 
Efficiencies for massive subhaloes can in general be extracted using 
low resolution simulations. On the other hand, efficiencies for low 
mass subhaloes can be derived using a high-resolution simulation 
of low mass objects. 



6.3 Why do SAM & HYD agree? 

In the previous sections we have shown that our SAM agrees quite 
well with the HYD, even when comparing individual galaxies. This 
fact seems to be surprising. After all, the HYD follows the dynam- 
ics of typically 10 3 — 10 6 particles within a galaxy, and the rates 
of cooling, SF, and feedback should be affected by many details 
within the HYD. For example, metallicity, the abundance of dif- 
ferent elements within the gaseous halo, and t he gas density p ro- 
file should all affect the cooling rate strongly jMo et al.ll2010l) . It 
is also known that the morphology of the galaxy, and the spatial 
distribution of cold gas w ithin the disk, should affect the SF rates 
(iMcKee & O striker! 12007). How is it that the simple SAM, based 
on average efficiencies per subhalo mass and time, can reproduce 
all these dependencies accurately? 

The fact that the simple, one-phase model, provides an accu- 
rate match to individual galaxies within the HYD might help us to 
discover the reasons for the good agreement. We will thus build 
our arguments using the one-phase model, since it is simpler, and 
includes an analytic integral solution. 

Assume that for a specific galaxy within the HYD, the values 
of f s deviate from the average values as adopted by the SAM. From 
Fig. [12] it looks like deviations between different merger-trees are 
negligible. This means that within the full history of one galaxy 
(including all its progenitors), the deviations in / s with respect to 
the SAM efficiencies should average to zero. In order to understand 
how these deviations affect ra s tar we examine the analytic solution 
presented in Eqs.[l9]&|20] As seen from these equations, the mass 
components of the galaxy depend on integral values, such as the 
integral of / s over time. Since the deviations in / s should be aver- 
aged out over time, these deviations do not affect ra sta r much. This 
is also true if the model is not ideal in the sense that / s depends on 

The fact that our efficiencies are computed using average val- 
ues per subhalo mass and time is the other key ingredient for the 
good match between the SAM and the HYD. Such a method guar- 
antees that on average, the mass of each baryonic component within 
the SAM will agree with the HYD. Moreover, the scatter between 
the models should be small because various random processes can- 
cel each other: the SF history of one galaxy is built from different 
episodes in its history, and the stellar mass is a sum over all these 
episodes. If each SF event within the HYD includes some random 
error (with respect to the SAM), the sum of all should have a much 
smaller scatter. This is the reason why the deviations in SF rates 
between the HYD and SAM are of the order of 0.5 dex, while the 
total stellar masses agree much better, to a level of 0.1 dex. 
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7 SUMMARY AND DISCUSSION 

In this work we have developed a method to post-process a hydro- 
dynamical simulation of galaxy formation (HYD), and to extract 
the simple baryonic laws that shape the evolution of galaxies within 
it. By using the same laws within a semi-analytic model (SAM), we 
confirmed that the resulting galaxies in both models remain almost 
unchanged. 

We have used a state-o f-the-art HYD, taken from the OWLS 
project dSchave et all l201Ch . This simulation includes radiative 
cooling, a galactic wind model for feedback, and a SF recipe that 
mimics the observed Kennicutt- Schmidt law. On the other hand, 
we have used a simple SAM, based on the approach presented in 
Neistein & Weinmann ( 2010). In this SAM, the processes of ac- 
cretion, cooling, SF, and feedback, are modeled using efficiencies 
that depend only on the host subhalo mass and redshift. Although 
previous studies did not find a good agreement between HYDs and 
SAMs, our simple SAM can reproduce the results of the advanced 
HYD at the level of 0.1 dex, when comparing individual galaxies 
(see Fig. [7]). The same level of accuracy is achieved for all red- 
shifts and for various mass components (stars and gas, although 
the cold gas component shows deviations of ~ 0.2 dex). Statistical 
properties of galaxies, like the stellar mass function, are accurately 
reproduced by the SAM as well (see Fig.fTOt. 

We claimed that the main reason for the agreement between 
the SAM and the HYD is the nature of the differential equations 
that govern galaxy evolution: Using a simple one-phase model, 
that can be solved analytically, we show that the mass compo- 
nents depend on the integrals of efficiencies over time, so temporal 
variations cannot be seen in the total stellar mass of galaxies. On 
the other hand, instantaneous properties like SF rates, reveal much 
larger deviations between the SAM and the HYD, at the level of 
0.5 dex. The fact that our recipes are matched for the same subhalo 
mass, implies that objects from the SAM and HYD should agree on 
average, since the common reference between these models is the 
merger-tree, based on the subhalo mass. 

We have found that satellite galaxies within the HYD experi- 
ence stripping of hot gas, on a time scale of ^ 3 Gyr. Although we 
have used only one time- scale for stripping, we found that a more 
accurate description would demand a time- scale that depends on 
the satellite mass, and on the infall redshift. In addition, satellite 
galaxies have different efficiencies of cooling and feedback than 
central galaxies, a fact that was neglected in this work. Due to the 
above simplifications made in our model, the agreement between 
satellite galaxies in the HYD and SAM reaches a value of 0.2 dex 
at low redshift. 

The method presented here is very robust and should work 
for any hydrodynamical simulation. We have shown that it even 
works if we assume only one phase of gas within subhaloes. In 
this highly simplified SAM there are only two processes: smooth 
accretion, and SF. Nonetheless, the SAM is able to reproduce the 
same galaxies as in the HYD with the same accuracy as a SAM 
with more phases of gas (see Fig. [13]) • This test proves that our 
method is insensitive to the specific way the SAM equations are 
written, and that there are many possible models to describe the 
same simulation. 

We have shown that the variance in the baryonic efficiencies is 
very small, once efficiencies are first averaged over the progenitors 
of each z = galaxy. We have also extracted the SAM efficien- 
cies from the progenitors of one random massive galaxy at z = 0. 
These efficiencies are then able to reproduce the full population 
of galaxies within the HYD. Consequently, our method can sepa- 



rate between baryonic processes and cosmic variance, and can be 
used to interpret one zoomed hydrodynamical simulation. In a fu- 
ture work, we plan to use the method developed here in order to 
scan combinations of various implementations of feedback and SF 
within hydrodynamical simulations. 

We have investigated the efficiencies of the baryonic processes 
within the HYD, showing that this specific simulation deviates 
from the values adopted by standard SAMs. We have shown that 
smooth accretion does not always fol low dark-matter, and is le ss ef- 
ficient for low-mass haloes (see also Ivan de Voortetai1l201ll) . All 
the fresh gas falling into subhaloes is either hot or dilute, and can- 
not form stars. For a given subhalo mass, cooling time-scales are 
roughly proportional to the cosmic time, with somewhat shorter 
time-scales at high-z. On the other hand, the efficiency of SF does 
not change significantly with cosmic time and does not decrease 
strongly for high-mass haloes. Lastly, the HYD used here does not 
show noticeable contributions from merger-induced star formation 
bursts. We hope that these findings can help to bridge the gap be- 
tween SAMs and HYDs, and will be useful for interpreting various 
existing models. 
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APPENDIX A: MORE EFFICIENCIES 

In this section we plot additional efficiency values that were ex- 
tracted from the HYD. In Fig. IA1I we plot the smooth accretion, 
when stripping is allowed for central galaxies. These accretion rates 
are very different from our standard model (Fig. O, proving that 
stripping occurs often in low mass haloes. 

In Fig. IA2I we show the cooling and feedback efficiencies 
when we use an additional gas phase, ra e j c to describe the ejected 
gas. Here f c describes the transition from rahot to ra co id only, and 
is normalized by rahot- Feedback is defined as the total heated 
gas, including all gas particles that started within ra co id and ended 
within rahot or ra e j c . 

The efficiencies plotted in Figs. [3ll6l use a fixed particle mass, 
and do not account for mass loss due to SN and stellar winds. Al- 
though these efficiencies were treated self-consistently in the main 
body of this work, they might deviate from the true efficiencies, 
computed with the proper particle mass. In Figs. lA3l & IA4l we plot 
the true efficiencies. These are very similar to the values used be- 
fore, but show some small differences in the overall normalization. 
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Figure A2. The cooling and feedback efficiencies within the HYD, using 
a separate ejection phase. The definitions of line-types are the same as in 
Fig.0 
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Figure A3. The cooling and feedback efficiencies using the proper mass for 
each particle within the HYD. The definitions of line-types are the same as 
in Fig. [3] 



Figure A4. The accretion and SF efficiencies using the proper mass for 
each particle within the HYD. The definitions of line-types are the same as 
in Fig. [3] 
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